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ABSIHACT 


An energy audit carried out by the Tata Energy Research Institute shows muffle furnace used 
in the Firozabad glass industry operate at very low efficiencies, of the order of 2 %. 
Considerable savings in fuel costs would be achieved if the energy efficiency of these muffle 
furnaces is increased even moderately. Hence the objective of this research has been to develop 
a mathematical model which can be used for the energy-efficient design of the Firozabad muffle 
furnaces. The mathematical model involves the simultaneous solution of Bernoulli equation and 
energy balance equations. Bernoulli’s equation has been used to determine the mass flowrate of 
gas inside the muffle furnace by natural draught. Computations showed that the calculated hearth 
and flue gas tempearture did not match well with measurements. The discrepancy between 
predicted and measured temperatures can be attributed to the low mass flowrate of flue gas 
within the furnace which implies that the assumption of constant hearth temperature in our model 
may not be appropriate. 
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NOMENCLATURE 


m 

C P 

Cf 

V 

h„ 


l coal 


: hl 

T h 2 

T h3 

T 0’ T N1 

T N1’ T N2 

T N2 ,T N3 

T N3 ,T N4 

T N4 ,T N5 

T N5’ T N6 

T N6 ,T N7 

T N7 ,T N8 


Q 

D 

NX 

NY 

NC 


mass flow rate 
specific heat 
correction factor 
volume of coal bed 

convective heat transfer coefficient across the vertical walls. 

thermal conductivity 

temperature of coal bed 

ambient temperature 

temperature inside the bottom hearth 

temperature inside the middle hearth 

temperature inside the topmost hearth 

temperature at the extremeties of flue gas stream 0, 1 

temperature at the extremeties of flue gas stream 1,2 

temperature at the extremeties of flue gas stream 2,3 

temperature at the extremeties of flue gas stream 3,4 

temperature at the extremeties of flue gas stream 4,5 

temperature at the extremeties of flue gas stream 5,6 

temperature at the extremeties of flue gas stream 6,7 

temperature at the extremeties of flue gas stream 7,8 

cross-sectional area 

Heat lost or heat gained 

Thickness 

Number of divisions in the horizontal stream (0-1, 2-3, 4-5, 6-7) 
Number of divisions in the vertical stream (1-2, 3-4, 5-6) 

Number of divisions in the vertical stream flowing through the chimney 



SUBSCRIPTS 


01 

related to flue gas stream 0,1 

12 

related to flue gas stream 1 ,2 

23 

related to flue gas stream 2,3 

34 

related to flue gas stream 3,4 

45 

related to flue gas stream 4,5 

56 

related to flue gas stream 5,6 

78 

related to flue gas stream 7,8 

chm = 

related to chimney 

coal = 

related to coal 

hi 

related to bottommost hearth 

h2 

related to middle hearth 

h3 

related to topmost hearth 



1. INTRODUCTION 


A furnace is a device in which the chemical energy or the electrical energy is converted into heat 
which is used to raise the temperature of materials. 

The function of an industrial furnace is to generate and apply heat in such a manner that the 
heated product will conform to certain specifications, with lowest cost of heating per unit of finished 
products. 

The selection of that source of heat which is best suited to a particular case necessitates the 
possession of the knowledge of - important properties of fuel, the equipments needed for their 
preparation the furnace for burning fuel economically with the correct amount of air, which varies with 
the nature of fuel, the devices for controlling furnace temperature and furnace atmosphere, cost of a 
given amount of heat when produced by different fuels. 

The "economy" or "efficiency" when used in their true sense in connection with industrial 
furnaces, have reference to heating cost per unit weight of finished saleable product. "Heating cost" 
includes not only the cost of fuel but also the cost of firing the furnace, the cost of maintainence and 
repair, the cost of generating a protective atmosphere and cost of burnt, spoilt or rejected pieces. 

With so many different items entering into the cost of heating, it is quite possible that in some 
cases, the highest priced fuel may be the cheapest in end, so far the total cost is concerned. With 
favourable conditions, excellent design, and with good operation appreciable furnace efficiencies can 
be obtained. By efficiency, it is meant the fuel efficiency which is the ratio of heat input into the stock 
to the potential heat of the fuel. Gases, for example, can give up heat to the charge only as long as they 
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are hotter than the charge. Consequently, the flue gases leave industrial furnace at a very high 
temperature. This factor is responsible for low furnace efficiency. 

It is this requirement of low heat cost per unit of heated products that has led to the development 
of concept of the energy audits. These audits; take into account the various factors, (discussed in 
background) leading to heat loss and hence affecting the economy of the process adversely. One such 
audit carried for a typical muffle furnace at the Firozabad glass industry, by the Tata Energy Research 


Institute reveals the following energy balance: 

Flue gas losses 

= 47.8% 

Structural losses 

= 36.0% 

Heat loss due to CO formation 

= 9.9% 

Heat loss to trays 

= 3.7% 

Heat loss to bottom ash 

= 1.1% 

Heat gain by bangles 

= 1.5% 


The above result implies that the efficiency of the muffle furnace is as low as 2%. The existence 
of the scope for a considerable improvement, via the design of an energy efficient muffle furnace, and 
hence deriving the maximum economic benefits, has inspired the project undertaken. 
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2. BACKGROUND 


2.1 HEAT DISTRIBUTION IN A FURNACE 

Thermal efficiency of a furnace which is the ratio of the quantity of heat utilised to the quantity 
of heat supplied, varies in a wide range from as low as of 10% in some furnaces to 60% or even more 
than 90% in case of a well designed one [1]. 

The fuel required for combustion is burnt in the combustion zone of the furnace. Contrary to 
the requirements, not all the heat is absorbed by the stock. Some of it passes into the furnace wall and 
hearth. Another portion of heat is lost to the surrounding by radiation and convection from the outer 
surface of the walls or by conduction, into the ground. Heat is also radiated through cracks or other 
openings and furnace gases pass out around the door, frequently burning in the open and carrying off 
the heat. Heat is also lost each time the door is opened. Finally, the bulk of heat loss passes out along 
with the flue gases either in the form of sensible heat or as incomplete combustion. 

Various factors which affect the fuel economy and furnace efficiency adversely are: 

Improper design 
Incomplete combustion 
Improper heat distribution 

Operation at a temperature much higher than desired temperature. 

Heat losses from walls, openings, cracks. 

Improper control of furnace draught 

No waste heat recovery from the furnace flue gases. 
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The present work is primarily concerned with heat balances. Hence the various kinds of heat 

losses ?ire considered. Heat losses can be classified into following four categories: 

1. Heat loss from furnace walls - Heat is lost to the surroundings by radiation and convection from 
outer surface of furnace walls and by conduction into the ground. 

2. Heat loss through openings and crack - There may be various openings in the furnace, for 
example,peep holes, inspection doors. Charging and discharging doors, which cause considerable 
heat loss. 

3. Heat loss through waste gases - In combustion type furnaces two additional sources of heat loss 
occur. These losses are caused by the heat energy that the product of combustion take out of 
the furnace in the shape of unburnt fuel and in the form of sensible heat due to high temperature. 
The heat loss is commonly known as stack loss. 

4. Special losses 

a) Heat losses due to part of stock projecting out of the furnace - Heat is lost due to 
dissipation of heat from metal portion projecting from the inside of the furnace into the 
air. Throughout the heating period, heat flows from the part that is in the furnace to the 
part that is outside, and a portion of this heat is dissipated from the latter part, by 
radiation and convection. 

b) Heat loss to tongs and charging machines - Cold tongs absorb heat while in the furnace. 
Although the flow of heat into tongs is less than the heat that is radiated through the open 
door, it is by no means inconsiderable. 

c) Losses to trays, conveyer chains, and rollers. 
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2.2 DETAILS OF FIROZABAD GLASS FACTORY 


The Firozabad glass factory accounts for nearly seventy percent of the glass production in the 
small sector in India. The cluster is an agglomeration of small units engaged in the manufacture of 
hollow ware, headlight covers, beads, decorative items and bangles. Principally, three types of furnace 
are employed by the Firozabad industry - tank furnaces, pot furnaces, muffle furnace. The muffle 
furnace commonly known as ‘pakai bhatti’ is used for baking bangles to give them the required finish. 

There are about 400 registered muffle furnace units in the cluster. These furnaces work round 
the clock and for nearly 350 days a year. Coal is the primary fuel for these furnaces although wood 
is also used as a supplementary fuel in some of the units. Due to the irregular supply and sometimes 
the non-availability of coal from government quota, coal is generally purchased from open market at 
prices in the range of Rs. 2500-3000 per tonne. The schematic diagram of a typical muffle furnace used 
in Firozabad is shown in Fig. 1. Around 14,4000 bangles per day are baked in a single furnace. The 
bangles are baked in trays through three sequential baking hearth. The temperature in these are around 
400,605,785 degree Celsius in the top, middle, and bottom hearth, respectively. Bangles arranged on 
a steel tray are first baked on the top hearth for 55-60 seconds (low temperature); on middle hearth for 
55-70 seconds; and finally on the bottom hearth for 85-105 seconds (high temperature). There is no 
temperature measuring device and the actual baking time is decided by the operator based on his 
experience. 

The muffle furnace operates on the natural draught. Fuel (coal/wood) is fed over the fixed grate 
provided at the bottom of the furnace. Hot gases pass through a zig-zag path before leaving the 
chimney. The flue gases flow upwards around the hearth finally leaving the furnace at around 880 
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degree Celsius. Heat is transferred from the flue gases to the hearth by conduction. The present cost 
of construction of one muffle furnace is about rupees 7000/- [2], 
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3. PROBLEM FORMULATION 

3.1 ASSUMPTIONS 

The construction details of the furnace and the composition of the coal was not available. Hence 
the following assumptions have been made. 

1) The coal used for the purpose is lignite. The composition (% weight) of which is : 

C = 59.9% 

H 2 = 4.37% 

0 2 = 18.64% 

N 2 = 1.22% 

S = 2.65% 

Ash = 13.22% 

The adiabatic flame temperature is 2060 degree Celsius [3]. 

2) The actually obtained flame temperature is about 70% of the theoretical attainable 
temperature [4]. 

3) The coal particles are assumed to be perfect spheres with a diameter of about 40mm. 

4) The coal pieces are so arranged that it forms a close packed bed. 

5) The air enters the furnace at an ambient temperature (40 degree Celsius). 

6) The process is assumed to be in steady state. 
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3.3 HEAT BALANCE EQUATIONS 

The temperature should increase from T 0 to T N1 due to the net inflow of heat from combustion. 
Thereafter, it is expected to decrease because of losses to the hearth and the ambient. 

The energy balance of the various air streams results in the following equations: 

Stream 0-1 
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Implicit in the equation 1 is the assumption that gas is in contact with coal along the entire length of 
stream 0 , l . In the expression for Q coa |, an adjustable parameter, Cf, has been incorporated to account 
for the lack of information about (a) properties of the coal, (b) the nature of gas-solid contact on the 
grate and (c) heat losses to the ground and through furnace wall below the grate. Rj, R 2 , R 3 , R 4 
represents the resistances (per unit interfacial area) to heat flow from the flue gas to different parts of 
furnaces, as shown in Fig. 2. Implicit in equation 1 is the asumption that the inner walls of the bottom 
most hearth is at a uniform temperature, Tj_j j. A simple energy balance on a differential element in 
stream 0,1 will show that the flue gas temperature varies exponentially. Hence the linear variation in 
flue gas temperature assumed in equation 1 can only be valid if each gas stream is divided in a number 
of small element (NX). Equation l has to be written for each element in stream 0,1. A similar 
methodology has been followed for carrying out energy balance for all other air stream shown in 
Fig. 2 . 

Stream 1-2 
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Stream 3-4 
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where, 
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Three more equation can be obtained by doing a steady state heat balance on each of the three hearths 
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Hearth - 2 
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Equation of the form 1 to 1 1 are then solved simultaneously after dividing the stream into an 
"adequate" number of subsections. 
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4. CORRELATIONS 


The various correlations used for the calculation purposes discussed below. The solution of 
equation 1 to 1 1 requires a prior knowledge of the gas mass tlow rate. It can be obtained with the help 
of Bernoulli’s equation, as shown below. 

4.1) Bernoulli’s equation for the calculation of mass tlow rate [5): 


Ci) 

aU 



Fire Zone 

Fig. 5 Natural convection in a chimney. 


As air heats up in a fire place it expands so that the gas flowing inside the chimney stack is lighter than 
the air outside. The pressure at the lower and upper ends of the stack are determined by the relatively 
denser outside air, so that 

(P 2 -P 1 )= -pj gH 


17 


where. 


H = height of the chimney stack 

pj = density of the cooler air outside 
P 2 = density of the warmer air inside 

The lighter air inside the chimney stack, therefore experiences buoyancy and moves upwards. 
Since the density of fluid is not constant everywhere, the mechanical energy equation cannot be applied 
between regions 1 and 2. The problem is simplified, if the flow is divided into three regions. The first 
region extends from l to 3 where the colder air flows at constant density p^ from the relatively stagnant 
region into the fire zone of cross sectional area A^. The second region is the fire zone itself where heat 
is added to the air such that its temperature rises at essentially constant pressure. In the third region 
from section 4 to 2 the lighter gas flows (at constant density pj) into the stack of area A 2 finally exiting 
into the atmosphere at pressure P 2 . The mechanical energy equation applied to the first and third region 
to give (on neglecting losses): 


£1 

Pi 




and 



The continuity equation relates the velocities at sections 2, 3 and 4 as 
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Expressing all the velocities in terms of the desired stack velocity V 2 , also, P 3 = P 4 and 

P 2 -P x = -p lg H 


we get 


K - 




2(— “l)gH 

Po 


A -y , P, 
^1 Pi 


The mass flow rate is calculated as shown: 


m = P y4 2 V 2 


4,2) Calculation of heat transfer coefficient for the packed coal bed 

The coal is burnt below the bottommost hearth. As assumed, the perfect spherical coal particles form 
a packed bed. The heat transfer coefficient for the same is given by [ 6 ]: 

,03 

h = 12 *v x (- — ) 

°>r K .us' 

d p 
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where. 


h v = hs = volumetric heat transfer coefficient, kcal/cubic metre-hr degree Celsius. 

t = temperature in degree Celsius. 

dp = particle diameter, mm. 

vog = superficial gas velocity, m/sec 

4.3) Calculation Of Heat Transfer Coefficient across vertical plates 

For the free convection across the vertical plane surfaces, the correlation for the plate length 
Nusselt number is given by 

Nu, = CRa, m 

where the coefficient C and the exponent m are functions of Raj, the rayleigh number based on the total 
plate length. In laminar region (10‘*<Raj< 10 9 ), m = 1/4. For the turbulent region (10 9 <Raj) m = 

1/3. 

The correlation that has been used is [7]: 

Nuj = 0.68 + 0.670 (Raj 0 - 25 ) [1 + ((0.492/Pr) 9/ 1 6 )](" 4/9 > 

for 0<Raj< 10 9 

Nuj = (0.825 + 0.387 (Raj 1/6 ) [1 + ((0.492/Pr) 9/16 )]^ 8/27 )} 2 

for 10 9 <Raj 


20 



5. 


RESULT & DISCUSSION 


algorithm 

1) Assume values for CF, XN, ZN and CN 

2) Guess the value for T^g , temperature of the flue gas coming from the stack. 

3) Calculate mdot using Bernoulli’s equation. 

4) Calculate temperature profile in the furnace by solving equations 1 to 11 
for each element in the flue gas stream. 

5) If convergence in temperature profile not achieved then go to step (3) using 
the new value of Tj^g. 

Calculation is carried out for the furnace depicted in fig. 1 . The internal dimensions of the 
furnace such as (a) dimensions of the flow channel, (b) hearth dimensions and (c) wall thickness are not 
known. Some of these information have been estimated from fig. 1. by assuming that the furnace has 
been drawn to scale. 

Initially calculations were carried out by dividing each section of the flow stream into 10 
sections(XN=ZN=CN = 10).The ambient air ,and hence T is assumed to be 40 degree Celsius. The 
value of CF was varied until the value of Tjqg =850 degree Celsius was obtained. It may be recalled 
that this is the measured flue gas exit temperature. Table 1 shows the temperature distribution inside 
the furnace at CF=0.04 when TN8=850 degree Celsius is achieved. 
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Table 2: Temperature profile at CF = 0.04, XN = ZN = CN =50 


Sr. No. 

Symbol Used 

Temperature in U C 

1 

T 0 

5745.42 

2 

T N1 

3254.61 

3 

T N2 

2911.06 

4 

T N3 

2148.68 

5 

T N4 

2656.78 

6 

T N5 

2709.27 

7 

T N6 

2470.20 

8 

T N7 

1172.88 

9 

t hi 

3564.43 

10 

T H2 

2384.02 

11 

T H3 

3044.68 


MASS FLOW = 1.855 * 10' 3 
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Using XN = ZN= CN = 200, also does not significantly alter the nature of results shown in table 
1 and Table II clearly using Tj^g has a convergence criteria is not desirable because it lead to 
unusually high value temperature. The measured hearth temperature were 400, 605 , 785°C in the 
top, middle, and bottom hearth respectively. Moreover it is also evident that the predicted 
temperature profile does not decrease continuously as expected and let the temperature of the 
topmost hearth is greater then the metal hearth. The second approach for solving this problem is to 
use a value of CF such that the calculated hearth temperature are similar to the measured values. 
Table III and IV depict the temperature profile for CF = 0.007 when each stream is divided into 
10 and 50 section respectively. 
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Table 3 : 


Temperature profile at CF = 0.007, and XN = ZN = CN = 10 


Sr. No. 

Symbol Used 

Temperature in 

1 

T 0 

1374.80 

2 

T N1 

763.975 

3 

T N2 

603.265 

4 

T N3 

434.578 

5 

T N4 

488.373 

6 

T N5 

481.532 

7 

T N6 

440.968 

8 

T N7 

285.885 

9 

t hi 

770.594 

10 

T H2 

460.722 

11 

T H3 

535.754 


MASS FLOW = 2.895 * 10~ 3 
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Table 4: Temperature profile at CF = 0.007, XN = ZN = CN =50 


Sr. No. 

Symbol Used 

Temperature in U C 

1 

T 0 

1473.06 

2 

T N1 

878.159 

3 

T N2 

764.323 

4 

T N3 

598.351 

5 

T N4 

737.095 

6 

T N5 

763.069 

7 

T N6 

700.382 

8 

T N7 

432.864 

9 

t hi 

912.145 

10 

T H2 

648.375 

11 

T H3 

856.491 


MASS FLOW = 

= 2.726 * 10' 3 
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Even though the temperature profile of the flow stream for XN = ZN = CN = 50 seems to be 
more realistic, the middle and top hearth temperature are inconsistent with observations. 

To get realistic value of temperatures, the procedure for solving the hearth temperature was 
modified. Instead of solving the three hearth temperatures T H1 , T H2 & T H3 , using the energy 
valance approach as manifested in equations 9-11, to TH3 was calculated as follows : 


T H i- 0.1 ( T + t n1 + t N2 + T N3 ) 


T H2 ~ °- 1 (T N2 + T N3 + T N4 + T N5) 


T H3 ~ °- 1 ( T N4 + t N5 + T N6 + T N7) 


A coefficient of 0. 1 was used in the equation instead of 0. 1 to ensure the that the hearth temperature 
was lower than that of the gas stream flowing around it. Unfortunately, this scheme also did not 
provide the desired results. 

The main reason for improper prediction of temperature seems to be the low mass flow rate of flue 
gas inside the furnace. Lower mass flow rates imply greater heat dissipation to the surrounding and 
hence a greater drop in temperature of the flow stream, a fact that can be easily verified the energy 
balance equations listed in the previous chapter. One method of overcoming w in problem is to 
divide the stream into greater number of division. Unfortunately, dividing each stream into as many 
as 200 section, that is dividing the entire flow path into the furnace into 1600 section did not 
rectified the situation clearly the problem lies elsewhere. 


27 



It is seen that the most plausible cause for the factors to predict the temperature profile in the muffle 
furnace is the assumption that all the hearth walls are at a uniform temperature. At low mass flow 
rate the temperature of the flue gas would be expected to be in equilibrium with the hearth walls . 
That is the temperature of the hearth wall should very enough manner similar to the flue gas stream 
adjacent to it, although its temperature would be lower than the flue gas temperature. For example 
the temperature of the hearth will just above point ) would be the highest and it would decrease 
continuously upto N3. A similar situation should exist in the other two hearths. 

Therefore assumption of uniform hot temperature resulted in larger heat loses from the flow stream. 
Hence we have a situation where heat flows from, example a hearth 1 to stream 2 - 3 (b) hearth 2 to 
stream 4-5 and (C) hearth 3 to stream 7-8. 

Future Work 

The energy balance equations stated in the previous chapter has to be modified to account for the 
non-uniform hearth temperature. 

(1) Equation 1 to 8 should be modified & make temperature Tjqj to as variables 

(2) Replacing equation 9 to 11 with "radiation’s" equation. 
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CONCLUSIONS 


The main reason for improper prediction of temperature seems to be the low 
mass flow rate of flue gas inside the furnace. Lower mass flow rates imply greater 
heat dissipation to the surrounding and hence a greater drop in temperature of the 
flow stream , a fact that can be easily verified by the energy balance equations listed 
in the previous chapter. One method of overcoming this problem is to divide the stream 
into greater number of divisions. Unfortunately, dividing each stream into as many as 
200 sections, that is dividing the entire flow path in the furnace into 1600 section did 
not rectify the situation. Clearly the problem lies elsewhere. 

It is seen that most plausible cause for the factors to predict the 
temperature profile in the muffle furnace is the assumption that all the hearth walls 
are at a uniform temperature. At low mass flow rate the temperature of the flue gas 
would be expected to be in equilibrium with the hearth walls. That is the temperature 
of the hearth wall should be vary in a manner similar to the flue gas temperature. For 
example the temperature of the hearth will just above point 0 would be the highest 
and it would decrease continuously upto N3. A similar situation should exist in the 
other two hearths. 

Therefor assuption of uniform hearth temperature resulted in larger heat losses from 
the flow stream. Hence we have a situation where heat flows from, example a hearth 
1 to stream 2-3 (b). hearth 2 to stream 4-5 and (c) hearth 3 to stream 7-8. 



FUTURE WORK 


The energy balance equations stated in the previous chapter has to be modified to 
account for the non-uniform hearth temperature. 

(1) . Equation 1 to 8 should be modified & make temperature Thi to Th 3 as variables. 

(2) . Replacing equation 9 to 1 1 with “radiation’s” equation. 

% 

(3) . The problem should be worked under steady state conditions. 

(4) . The temperature change should be two-dimensional. So differential equations 
have to be solved. 



REFERENCES 




[1] Gupta, O.P., Elements of fuels, furnaces and refractories, Khanna Publishers, 1990. 

t 

[2] Pamphlet from Tata Energy Research Institute 

[3] Trinks, W., Industrial Furnaces, Volume II, John Wiley and Sons, Inc., Chapman & Hall, 

Ltd., London. 19 73 

[4] Schuhmann, R., Metallurgical Engineering, Addison - Wesley Publishing Company, Inc. I9&I 

[5] Gupta, Vijay and Gupta S. K., Fluid Mechanics & Its Applications, Wiley Eastern Ltd J 992 

[6] Gieger, G.II. and Fourier, D.R., Transport Phenomenon in Metallurgy, Addison- 
Wesely Publishing Company. J972 

[7] Chapman Alan, J., Heat Transfer, Macmillan Publishing Company, New York, Collier 
Macmillan Publisher London. 15 9 9 


29 



